A1 = readmatrix('NiI2_MT_COOL.dat')

x1 = A1(:,3)
y1 = A1(:,59)/0.003

hold on
%scatter(x1,y1,5,[0 0 0.8],'filled')
plot(x1,y1,Color=[0 0 0.8])
ylabel('Magnetization [emu/g]')
ylim([0.015,0.021])
xlim([40,80])
xticks([40:10:80])
xticklabels('')

set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off
ax = gca;
ax.XColor = 'white';
ax.YColor = 'white';

set(gcf,'Position',[0, 0, 292, 150])